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Monochromatic 7-rays are thought to be the smoking gun signal for identifying the dark matter 
annihilation. However, the flux of monochromatic 7-rays is usually suppressed by the virtual quan- 
tum effects since dark matter should be neutral and does not couple with 7-rays directly. In the 
work we study the detection strategy of the monochromatic 7-rays in a future space-based detector. 
The monochromatic 7-ray flux is calculated by assuming supersymmetric neutralino as a typical 
dark matter candidate. We discuss both the detection focusing on the Galactic center and in a scan 
mode which detects 7-rays from the whole Galactic halo are compared. The detector performance 
for the purpose of monochromatic 7-rays detection, with different energy and angular resolution, 
field of view, background rejection efficiencies, is carefully studied with both analytical and fast 
Monte-Carlo method. 

PACS numbers: 95.35.+d, 98.35.-a, 98.35.Gi 



I. INTRODUCTION 

The existence of dark matter (DM) in the Universe 
is widely accepted nowadays. The evidences come from 
many astronomical observations, which observed the 
gravitational effects of dark matter in different spatial 
scales, from dwarf galaxies, galaxies, galaxy clusters to 
the cosmological scale. It is recognized that DM particle 
should be neutral, cold and non-baryonic, which can only 
exist in theories beyond the standard model of particle 
physics. Among the large amount of DM candidates pro- 
posed in the literature, the weakly interacting massive 
particles (WIMP) are the most favored one, which can 
account for the observed DM density naturally. 

In order to determine the nature of the WIMP DM 
particles, we generally have three ways to probe the in- 
teraction between DM particles and the standard parti- 
cles: the direct detection measures the the scattering by 
DM with the detector nuclei; or produce DM particles di- 
rectly in a powerful collider such as LHC or ILC; finally 
the indirect detection searches for the DM annihilation or 
decay products in cosmic rays (CRs), including gamma- 
rays, electrons, positrons, protons, antiprotons and neu- 
trinos. For indirect detection the gamma-rays are usually 
the best probe of DM because they are not deflected by 
the magnetic field during the propagation. Further the 
technology of high energy gamma ray detector was devel- 
oped very fast in the last years. In general the gamma- 
rays are produced through the hadronization and decay 
of the DM annihilation/decay final states, and have con- 
tinuous energy spectra. Since there are large diffuse back- 
grounds of gamma-rays, it is usually not easy to figure 
out the DM signals from the background. Another anni- 
hilation channel is to monochromatic gamma-rays with 
small branching ratio via loop diagram. The monochro- 
matic gamma-rays by DM annihilation are usually taken 
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as the smoking-gun of the DM signal since there are no 
such background from astrophysical processes. If the de- 
tector has very good energy resolution, the background 
will be suppressed and the detectability will be improved. 

In this work we discuss DM annihilation into 
monochromatic photons, X0X0 — > 77 and X0X0 — ► 7^0 > 
with energy of the gamma-ray photon m x and m x — 
m 2 z /4m x respectively. Here we neglect the kinetic energy 
of DM particles since its movement is non-relativistic to- 
day. For very massive DM m x 3> mz, the photon energy 
of the two channels are identical and can not be dis- 
tinguished in experiments. This work tries to give the 
perspective of detecting such line spectrum gamma-rays 
from DM annihilation in the Milky Way, and show the 
requirements for detector design. 

The 7-ray flux from DM annihilation is proportional to 
the annihilation cross section and DM density square. As 
a typical WIMP DM we consider the lightest neutralino 
in the minimal supersymmetric standard model (MSSM) 
as an explicit example in our calculation. The cross sec- 
tion can be computed given the MSSM model parame- 
ters. In this work we will employ DarkSUSY package [1] 
to scan the MSSM parameter space. As to the Galactic 
DM density profile, numerical simulations indicate that 
DM is highly concentrated in the halo center. There- 
fore the Galactic center (GC) is usually the first choice 
searching for DM signals. In addition, there are also a 
large amount of substructures existing in the halo, mostly 
in the outer part of the halo. The contribution from 
substructures will also be discussed. The backgrounds 
for the monochromatic photon detection include CR nu- 
clei (mainly protons and Helium for energies < 10 TeV), 
electrons, diffuse continuous 7-rays and the 7-ray point 
sources in the GC region. The nuclei and electrons can 
be rejected through the particle identification technique 
of the detector. For the diffuse 7-ray background we need 
a high energy resolution to suppress the background. 

The paper is organized as follow. The MSSM model 
and the Galactic DM density distribution will be pre- 
sented in Sec. [TTJ In Sec. Mil the possible backgrounds 
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in our study are introduced, especially those in the GC 
region. Sec. IIVI gives the results of detection sensitiv- 
ity for the GC region through simple analytical estimate. 
We show in Sec. [V] the sensitivity for all-sky observa- 
tion in the scan mode, with substructures included, using 
Monte-Carlo simulation. Finally we draw the conclusion 
and discussion in Sec. IVII 



II. DARK MATTER ANNIHILATION INTO 
MONOCHROMATIC GAMMA-RAYS 



The 7-ray flux from DM annihilation can be written 



MS 



4>M = 
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(1) 



where i/j is a specified direction away from the Galactic 
center, p Q ps 0.4 GeV cm" 3 @, [3] is the local DM den- 
sity, Rq ps 8.5 kpc is the distance from the Earth to the 



GC, m x is the mass of DM particle, 



is the velocity 



weighted thermal average annihilation cross section, the 
multiplicity N = 1, 2 for jZq and 77 channels respec- 
tively. Finally J(il>) is the line-of-sight integral of the 
density square J( - 0) = p -i ^ / p 2 {l)dl. The 7-ray flux 



depends on both the particle parameters and the density 
distribution of DM. 



A. MSSM dark matter model 

Without losing the generality we take neutralino in 
MSSM as a typical WIMP DM in this work. In order to 
reduce the number of free parameters in MSSM we only 
take a few relevant parameters to our discussion free, as 
done in Ref. [J, that is, 



fx, M 2 , Mi, tan/3, M A , m , A b , A t , 



(2) 



where fi is the Higgsino mass parameter, M 2 and M\ are 
the wino and bino mass parameters respectively, tan /3 is 
the ratio of the vacuum expectations of the two Higgs 
fields, Ma is the mass of pseudo-scalar Higgs boson, too 
is the universal sfermion mass, A b and A t are the trilinear 
soft breaking parameters and the corresponding param- 
eters for the first two generations are assumed zero. 

We employ DarkSUSY to explore the parameter 
space of the phenomenological MSSM model |l]. The 
scan ranges of these parameters are: 50 GeV < 
\fx\,M2,Mi,M A ,m < 10 TeV, 1.1 < tan/3 < 55, 
sign(/i) = ±1, — 3too < A tl A b < 3mo. There are 
also other constraints from the theoretical consistency re- 
quirements and the accelerator data. Finally, we require 
the relic density of DM to be £l x h? < 0.128 according to 
the 3<T upper limits of WM AP seven year results [4| . 



B. Density distribution 

The most precise knowledge of the density profile of 
DM inside the halo comes from numerical simulations. 
Navarro et al. (1997) found that the density profile is 
universal for halos of different scales, referred as Navarro- 
Frenk- White (NFW) profile [| 



p(r) = 



(r/r s ) (1 + r/r s y 



(3) 



where p s and r s are two scale parameters depending on 
the mass and concentration of the halo. However, due to 
the limit of resolution, some other density profiles with 
different central behavior were also proposed in litera- 
ture. For example Moore et al. (1999) proposed the 
density profile with a much steeper inner slope [|| 



p(r) 
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The recent simulations tended to favor the Einasto pro- 
file with a gradual flattening of the logarithm slope of 
the inner behavior |7j. There were also studies showing 
that the density profile might be non- universal [8|. Con- 
sidering the diversity of the inner profile of DM density 
distribution, we adopt NFW and Moore profiles for this 
study. The model parameters are r s — 20 kpc, p s = 0.35 
GeV cm" 3 for NFW, and r s = 28 kpc, p s = 0.078 GeV 
cm" 3 for Moore profile respectively. The local density 
for these parameter sets is 0.4 GeV cm" 3 . 

To avoid the divergence of density when r — >■ 0, a cutoff 
scale is applied considering the fact that there should be 
a balance between the gravitational infall and the anni- 
hilation [9] . For common parameters of DM particle the 
maximum density is estimated to be ~ 10 18 M Q kpc" 3 



C. Substructures 

The cosmological structures form hierarchically in the 
cold dark matter scenario, that is, the DM collapses to 
form small halos first, then grows to larger and larger ha- 
los through accretion and merger. Numerical simulations 
show that there are a large number of subhalos surviving 
the merger history and existing in the Milky Way dark 
matter halo |ll| - [l4| . The dwarf galaxies are a part of the 
Galactic subhalos which have been observed. 

According to the simulations, the number density dis- 
tribution of subhalos as a function of its mass and loca- 
tion can be parameterized as [H, [lj| 



dn 



= v I 

dm suh ■ Anr 2 dr l + (r/r h ) 2 \ TO host 



(5) 



wit h ry, ps 0.14rhost for galaxy scale halo [16( and a ps 
1.9 (la . [l?]]. The normalization n is fixed by setting 
N{> 1O 8 M ) ps 100 (see [13 and references therein). 
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The density profile of each subhalo is also assumed to 
be NFW or Moore profile. To determine the density 
parameters of each subhalo, we use the concentration- 
mass relation given in Ref. [18| . The procedure of de- 
termining the parameters is as follows. For a halo with 
mass m su b, the concentration is derived according to the 
Csub - ""isub relation. Then we have r^ FW = r su b/c su b and 
r Moore r su b/0.63c su b following the definition of concen- 
tration [l8[. Finally p s is determined by the subhalo 
mass. 

To calculate the annihilation flux of photons from the 
subhalo population, we define the average density square 
as 



(Psub)O) = 



- j dm dn x f 

J su dm suh ■ 4:irr 2 dr J Vs 



2 

Psub 

V 5U b 



,dV, (6) 



where dV integrates over the volume of the subhalo. Re- 
placing p 2 in J(ip) in Eq. (Q]) with (p 2 uh ) we calculate 
the photon flux from DM annihilation in the subhalos. 

The J(V>) factor as a function of the angle ip away 
from the GC direction is shown in Fig. [TJ It can be seen 
that for the smooth halo the annihilation flux is highly 
concentrated in the GC. The substructure contribution 
is nearly isotropic in all directions. At large angles away 
from the GC the substructure component may dominate 
the annihilation flux for the Moore profile. 
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Figure 1: Angular distribution of DM annihilation luminosity 
(J factor) in the Milky Way. The inner plot uses log scale for 
x axis to show the details in the most central region. 



III. DETECTABILITY ANALYSIS 



The CR nuclei and electrons/positrons can be rejected 
through the design of particle identification technique 
of the detector, however, there still are a few particles 
which may be misidentified as photons and form the 
background. The combined nuclei fl ux, mostly proton 
and Helium nuclei, can be written as [19( 



ME) = 1.49 



E 



(~GeV 



-2.74 



cm 2 s 1 sr 1 GeV 



(7) 



which is an empirical formula from combined result of 
many measurements. For the electron plus positron spec- 
trum, we adopt a broken power law parameterization 



(t> e (E) = 1.5 x 1(T ] 
E 
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900 GeV 



cnr~ 2 s -1 sr _1 GeV _1 (8) 



according to the recent measurements by Fermi [20j and 
HESS [U In the following we employ two efficien- 
cies, r/ n and ry e , to represent the rejection power of the 
charged CRs. 

Then we come to the continuous 7-ray backgrounds. 
The first 7-ray background is the all-sky diffuse 7-ray 
emission, including Galactic and extra-galactic. For the 
extra-galactic 7-ray background we use the new measure- 
ment made by Fermi [23] 



(E) = 6.57x10 



7 (gcv) 



-2 - 

cm s 



^GeV 



( 9 ) 

The Fermi result of extra-galactic 7-ray emission is 
steeper than that obtained by EGRET [24]], which will 
result in an order of magnitude lower background when 
extrapolating to high energies (~TeV). 

Fermi collaboration also reported some data about the 
Galactic diffuse 7-ray emission (e.g., [23|, [2f| [26| ) , which 
is consistent with the results given by EGRET except 
the "GeV excess" [13]. Since the full Fermi data are un- 
available at the present stage, we use the EGRET data 
about the Galactic diffuse 7-ray emission in this work. 
The Galactic diffuse 7-ray flux is parameterized as [28j 



cj>f lac (E) =N (l,b)xlQ- 



where 



(cev) 



-2.7 



cm-V^r-^GeV -1 , 
(10) 



A. Backgrounds 

We first introduce the backgrounds which are essen- 
tial for the detectability analysis of the DM-induced 
monochromatic 7-ray signal. The backgrounds include 
charged CR particles, such as all kinds of nuclei, elec- 
trons and positrons, and continuous 7-rays. 



So. 5 



_ I V 1 + ( ( /35)V 1 + ( b /1.8) 2 



0.5 



85.5 



A /l + (;/35)V 1 + [ b /(l-l+0.022|i|)] 2 



|!| < 30° 
0.5 |Z| > 30° 



(11) 

in which the galactic longitude / and latitude b are ex- 
pressed in unit of degree. To make use of this measure- 
ment, we have to extrapolate Eq. (fTU|) to higher energies. 
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Besides the diffuse 7-ray emission, there are additional 
sources in the GC region. Since the GC region is very 
important for DM searches, the 7-ray sources in the GC 
region are necessary to be paid more attention. HESS 
observation showed there was diffuse 7-ray emission in 
the region \l\ < 0.8°J&| < 0.3° (GC ridge) on top of the 
diffuse background [29]. The spectrum is 



/GC-diff 

7 



(E) = 1.28x10 



- 4 f— V 

\GeVj 



cm s 1 sr 1 GeV 



. (12) 

Also there is at least one gamma-ray point source in the 
Galactic center area, which is labeled as 3EG J 1746-2851 
in EGRET catalog HESS J1745-290 in HESS catalog 
[U and OFGL J1746.0-2900 in Fermi catalog [11. The 
energy spectrum of this source given by HESS is [33] 



^point 



(E) = 2.5 x 10 



- 7 (— V 



cm- 2 s _1 GeV _1 , 



(13) 

which is valid between 200 GeV and 10 TeV. 

The differential fluxes of all the backgrounds men- 
tioned above are plotted in Fig. [2] Here we show the 
results within the circle with 1 degree radius around the 
GC. The electron and nuclei fluxes are multiplied by fac- 
tors 10~ 3 and 10~ 6 respectively, which represent the typ- 
ical rejection power of a 7-ray detector. We can see that 
the largest background is the diffuse source and the point 
source in the GC region. The Galactic diffuse 7-ray emis- 
sion is also important in the GC. However, for the sky 
regions far away from the GC, the extra-galactic back- 
ground would also become import. 
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all, we need particle discrimination to reject most of the 
charged CRs, say nuclei and electrons/positrons. This 
can be done through charge detection, neutron detec- 
tion and shower shape identification in calorimeter. The 
continuous 7-rays can be suppressed by improving the 
detector energy resolution. Furthermore, as can be seen 
in Fig. [1] the small region around the GC is the best 
candidate for the line spectrum 7-rays detection. Thus 
good angular resolution will also be effective to increase 
the signal-to-noise ratio. 

To simplify the study, we characterize the detector 
with some static parameters of the performance, that is, 
the parameters don't change with incident energy, di- 
rection, or particle type, like the energy resolution, the 
angular resolution, the field of view, the rejection power 
of electrons and nuclei and so on. The effective area of 
detector and exposure time are also key factors of the de- 
tection, however, the overall performance can be easily 
scaled given different effective area and exposure time. 
With the results of ideal detectors, we can provide the 
requirement of a real detector. 



C. Sensitivity calculation 



The event counts on a detector can be written as 



N i = r H -T eS -A-f i . 



(14) 



where rji is the detection efficiency of incident particle, 
usually assumed to be 90% for photon, and equal to the 
rejection power for electron and nuclei, T ef f is the effective 
exposure time for a given source or small sky region, A 
is the active area of detector, and /j is the flux of signal 
or background in given energy range and sky region. For 
a detector with energy resolution cr e , fi can be derived 
according to the differential flux <fi 

fi= d£l dE (t>i{E,l,b). (15) 

JAfl JE y ~3a E 

Here we choose the energy window to be ±3<te around 
E 7 , and AO, refers to the chosen observation field. Note 
for point source the above integral with respect to 57 dis- 
appears. 

The detection significance is defined as S — N s /^/T7b 
with N s the count of signal and Nf, the count of back- 
ground. The sensitivity can be derived as the minimal 
signal flux needed for a detection significance 5 = 5. 



Figure 2: Spectra of the possible backgrounds from the 1 
degree region around the GC. 



IV. GC REGION SENSITIVITY WITH 
ANALYTICAL ESTIMATE 



B. Ideal detector 

The sensitivity of monochromatic 7-ray detection is 
determined by the performance of a detector. First of 



In this section we estimate the sensitivity of monochro- 
matic 7-ray detection from the GC region by DM anni- 
hilation. We choose a typical planar calorimeter type 
detector with lm x lm size. With one year of flight on 
the orbit of International Space Station (ISS) and the 
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assumption of 90 degree field of view, the effective expo- 
sure time is ~ 7.9 x 10 6 s, which is about 1/4 of one year. 
More detail of effective exposure time can be found in 
Fig. E 

The background event number is 



Table I: GC detectability calculation configuration 
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(16) 



where AQ = Tr{od + I) 2 x (7r/180) 2 is the solid angle con- 
sidered for event direction selection, in order to keep the 
same number of signal events, where er^ is the detector 
resolution angle in unit of degree. Another method of 
events direction selection is to keep events of 1° from GC 
after point spreading of detector angle resolution. Con- 
sidering that the largest backgrounds come from point- 
like sources from GC, and they expands the same as the 
cuspy Moore dark matter profile, the one-degree selec- 
tion after expansion would result in same level of de- 
creasing of both signal and background events number, 
which means poorer sensitivity. On the other hand, if 
we select events with larger radius, we can keep basically 
the same number of signal events, and better sensitivity. 
Therefore the minimal monochromatic photon flux from 
DM annihilation within ~ 1 degree around the GC is 
/min = 5\fNb/rfr e sA, for a 5cr detection. 

The results of / m j n for different detector performance, 
i.e., different energy resolution, angular resolution, elec- 
tron rejection and nuclei rejection respectively, are shown 
in Fig. [3] The default detector performance settings are: 
energy resolution 1.5%, angular resolution 0.5°, electron 
rejection 10 -4 and nuclei rejection 10~ 7 . The different 
performances used in the calculation are summerized in 
Table U on page [5] The dots and triangles in the figure 
are the MSSM model predictions with a random scan in 
the eight-dimensional parameter space as introduced in 
Sec. Ill Al Both Moore and NFW profile are calculated 
and plotted. It is shown that the default detector config- 
uration is powerful enough to probe much of the MSSM 
parameter space for the model with DM density as cuspy 
as Moore type. If the DM density profile is NFW-like 
the sensitivity will be worse. We will discuss the effects 
of different density profiles more in the following. 

In Fig. [4] we show the comparison of sensitivities for 
longer exposure and larger area of the same type detec- 
tor. The upper limits of monochromatic photons from 
30 to 200 GeV derived according to 11 month Fermi data 
are shown by the line with errorbars [34| . Note the obser- 
vation region of Fermi's result is |6| > 10° plus 20° x 20° 
area around the GC, which is much larger than 1 degree 
around GC adopted here. Therefore the limits from GC 



common setup 


detector area: 1 m 2 , 
exposure time: 7.9 x 10 6 s , 
energy range: 
5nOeV~ lOTeV 


energy resolution 
angular resolution 
electron rejection 
proton rejection 


1%, 1.5%, 2%, 5%, 10% 
0.1°, 0.3°, 0.5°, 1.0°, 3.0° 
10°, 10 _1 , 10" 2 , 10" 3 , 10~ 4 
10~ 2 , 10" 3 , 10" 4 , 10" 5 , 10" 7 



would be smaller than the present given values. We just 
plot the results here for reference. 

In Fig. [S]the minimal detectable cross section is shown, 
for 77 (left) and jZq (right) channels respectively. In this 
plot the default detector performance is adopted, and we 
show the results for the NFW and Moore profiles, with 
different sky regions. It is shown that for the NFW profile 
the detectability will be orders of magnitude worse than 
that for the Moore profile. We also note that for the 
NFW profile the sensitivity is not affected significantly 
by the choice of sky area if the GC is included. How- 
ever, it is not the case for the Moore profile. This can 
be understood that for the Moore profile the signal flux 
drops very fast with the increase of the angle to the GC, 
which leads to a fast decrease of the signal-to-noise ratio. 
Therefore the large-angle average results in a worse sen- 
sitivity than the case focusing on the small region around 
the GC. The 95% exclusion limits of the cross section for 
the NFW profile from Fermi are also plotted (34[ . To 
compare with the capability of Fermi, we show the ex- 
clusion power of the detector with the default parameters 
mentioned before, at la level for the same density profile 
and observational sky region as Fermi. 



V. WHOLE SKY SCAN SENSITIVITY WITH 
MONTE-CARLO SIMULATION 

In this section we discuss the all-sky observation with 
scan mode. The effective exposure time is non-uniform 
at different directions for specific orbit of the detector. 
Therefore the simple analytical method discussed in the 
previous section does not work any more. We turn to use 
fast Monte-Carlo method to simulate the counts of signal 
and background, and calculate the sensitivity. 



A. Simulation configuration 

The flight orbit is assumed to be the orbit of ISS. We 
put an ideal detector on the orbit, and calculate the sig- 
nal and background particle counts by the Monte- Carlo 
sampling. We take the orbit data of ISS in 2003 as an 
example. The orbit of ISS has not changed too much 
in the period of several years. The sky-map of exposure 
time is plotted in Fig. [5] 



6 





Figure 3: The sensitivity with different detector features: energy resolution (top-left), angular resolution (top-right), electron 
rejection (bottom- left), nuclei rejection (bottom-right). The default settings are: energy resolution 1.5%, angular resolution 0.5°, 
electron rejection 10 -4 and nuclei rejection 10 -7 . The gray area at the bottom of each panel represents the flux corresponding 
to less than one event per year (exposure time) for a 1 m 2 detector. 



We choose detectors with a series of property param- 
eters, including energy resolution, field of view, electron 
and nuclei rejection powers. Seven neutralino masses, 
between 50 GeV and 1200 GeV, are chosen in the simula- 
tion. The details of the configurations are listed in Table 

IE 



Table II: Survey simulation configuration 



common configuration 


geometry factor: 3.14 m 2 sr, 
flight time: 1 year, 
flight orbit: ISS orbit 


simulated energy (GeV) 


77.6, 199.7, 416.3, 611.4, 
808.6, 1006.0, 1220.9 


energy resolution 
field of view 
electron rejection 
proton rejection 


1%, 1.5%, 2% 
40°, 60°, 90° 
10 -2 , 10" 3 , 10~ 4 
10" 5 , 10" 6 , 10" 7 



To get better statistics, a large enough cross section 
was chosen in the simulation. The photons from 77 and 
7Z0 annihilation channels are not distinguished, and only 
a sensitivity for single gamma line emission was given. 
Actually the photons from these two channels can not be 
discriminated when neutralino mass is larger than several 
hundred GeV, and if neutralino mass is small, there could 
be two peeks on the photon energy spectrum, thus the 
sensitivity would be higher. 

The overall contribution from DM annihilation, not 
only the smooth halo, but also the substructures, is taken 
into account. In the simulation we use Moore profile for 
both the smooth halo and subhalos. To compare with 
theory predictions, the model points are plotted along 
with the minimum detectable curves. In the calculation 
both Moore and NFW profile are used. And also the con- 
tribution of substructures and smooth halo are summed 
together. 
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with probability rj. 

The events kept in the efficiency selection are then re- 
constructed. The reconstruction of energy is simply im- 
plemented by a Gaussian smear of the inject photon en- 
ergy. The Gaussian width is er e , which represents the 
energy resolution of the detector. The nuclei rejection 
and electron rejection are simulated as a fixed pass rate 
as mentioned before. Passed electrons and nuclei are re- 
constructed the same way as photons. Figure [7] on page 
[5] shows a sample reconstructed energy spectrum for the 
neutralino mass 416 GeV. The direction reconstruction 
simulation is in the same way as energy reconstruction, 
with two-dimensional Gaussian smear. 



Results of simulation 



Figure 4: The sensitivity with different exposure factor and 
detector area. 



B. Event generation and reconstruction simulation 

The basic procedure of the simulation can be divided 
into two phases. First is the generation of the event sam- 
ples according to the angular and energy distributions of 
each component, including the signal and backgrounds. 
The second phase is to simulate the event reconstruct 
process based on the detector performances. 

For the signal events, the energy of generated photon is 
monochromatic for a given neutralino mass. The spatial 
distribution follows the result as shown in Fig. [1] The 
energy spectra and angular distributions of backgrounds 
are described in Sec. IIII Al In the simulation, we also 
need to know the direction of the detector. The detector 
pointing direction is set to be the same as the direction 
of ISS. 

The events are generated in interval of one second in 
detector operation time. To simplify the simulation code, 
static event rate for each component is used. The event 
rate is set to be the maximum value that can be re- 
ceived by the detector with acceptance Att ■ A, for energy 
range from 30 GeV to 1400 GeV, and a Poisson smear 
is used when event number in the one-second-interval is 
too small. However, the actual events received by a spe- 
cific detector will be much smaller than the above value 
since the effective area for photons with large incident 
angle will decrease. This effect is taken into account at 
the second phase through an efficiency factor r\. defined 
as ?/;Cos(p) when p < pfov, and when p > pfov, 
where p is the angle between the incident photon and 
the detector pointing direction, pfov is the maximum 
receiving angle of incident particles, and rji is the selec- 
tion efficiency for different incident particle. For photons 
the selection efficiency usually assumed to be 90%. For 
electrons and nuclei events, the detector rejection power 
is used as the selection efficiency. Then one event is kept 



The results of sensitivities for scan observation are 
shown in Fig. [5J for different detector performance. The 
area of detector is adopted to be 1 m x 1 m. The points 
represent the MSSM model predicted fluxes of the all-sky 
DM annihilation, with substructures included, for both 
Moore and NFW profile. Compared with the results of 
GC, the scan mode will be less sensitive for detecting 
the DM annihilation signal by more than one order of 
magnitude. 

VI. DISCUSSION AND CONCLUSION 

In this work we studied the strategy and feasibility of 
the monochromatic 7-ray detection from DM annihila- 
tion. We studied the detectability of the monochromatic 
7-rays from the GC and from the dark matter halo and 
subhalos, by analytical analysis and numerical Monte- 
Carlo simulation. The detector performance is studied 
for the purpose of monochromatic 7-ray detection. 

The sensitivity depends on the observation region and 
the DM density profiles. According to the Moore pro- 
file, the sensitivity of GC region would be much higher 
than other region. However, this is not true for NFW- 
like profiles. As we currently have little knowledge about 
the actual distribution, we can not rely on a particular 
model. Until now, high energy gamma observations, such 
HESS and Fermi, have not seen any evidence in the GC 
region. A safe strategy could be starting from the scan 
mode, then focusing on the interesting area if any evi- 
dence is found. On the other hand, the GC is always an 
interesting region for possible DM signals. The observa- 
tion of the GC should be tuned as much as possible. 

In order to detect the monochromatic 7-ray emission 
of DM annihilation in scan mode, the energy resolution 
and geometry factor, both detector area and field of view, 
are the most important detector feature. However, The 
resources are usually very limited on orbit, on both de- 
tector mass and power supply. With a given mass of the 
detector, thinner calorimeter would lead to larger active 
area. On the other hand, thinner calorimeters usually 





Galactic Longitude (°) 

Figure 6: Sky-map of effective exposure time according to 
the 2003 ISS orbit data. The whole sky area is divided into 
360 x 180 cells, with 1° x 1° of each cell. This map is derived as 
follows. The effective exposure time of each cell is calculated 
as cos(#) where 9 is the angle between the direction of the ISS 
and the direction of the cell. If 6 > 7r the effective exposure 
time is set to zero. Then the effective time of each cell is 
calculated in an interval of one second, and accumulated for 
all the time in one year (3.15 x 10 7 s). 

have poorer energy resolution. Therefore, to get better 
sensitivity, a careful balance is needed in detector de- 
signing. Also any detector design that would limit the 
field of view of the detector should better be avoided. At 
low energies (< 300GeV), the ability to reject cosmic-ray 
background, electrons and nuclei, are also important. To 
keep cosmic-ray background low enough, that is, lower 
than the major background, diffuse 7, at least 10 -3 elec- 
tron rejection and 10~ 6 proton rejection is required. 

If we are focusing on GC or other interesting region in 
pointing mode, the angular resolution, instead of field of 
view, becomes very important. In the GC region, the 7- 
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Figure 7: A sample reconstructed stacked photon energy spec- 
trum, for neutralino mass 416 GeV. 



ray emission is complex: several point sources exist there. 
Thus good angular resolution is needed to separate them 
apart, and suppress these point source backgrounds. Of 
course the energy resolution and detector area are still 
important features in order to detector monochromatic 
7-rays. Also, at low energy the rejection power require- 
ments of electron and proton are the same as in scan 
mode. 
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Figure 8: The minimum detectable flux of DM line emission in scan mode. Top- left: comparison among different maximum 
receiving angle (field of view); top-right: comparison among different energy resolution; bottom-left: comparison among different 
electron rejection power; bottom-left: comparison among different nuclei rejection power. The default parameters are energy 
resolution 1.5%, maximum receiving angle 90°, electron rejection 10 -4 , and proton rejection 10 _r . 
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